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Dengue  virus  has  traditionally  caused  substantial  morbidity  and  mortality  among  children  less  than  15  years  of 
age  in  Southeast  Asia.  Over  the  last  2  decades,  a  significant  increase  in  the  mean  age  of  cases  has  been  reported, 
and  a  once  pediatric  disease  now  causes  substantial  burden  among  the  adult  population.  An  age-stratified  sero¬ 
logical  study  (n=  1 ,736)  was  conducted  in  2010  among  schoolchildren  in  the  Mueang  Rayong  district  of  Thailand, 
where  a  similar  study  had  been  conducted  in  1980/1981 .  Serotype-specific  forces  of  infection  (7(f))  and  basic  re¬ 
productive  numbers  (R0)  of  dengue  were  estimated  for  the  periods  1969-1980  and  1993-2010.  Despite  a  signifi¬ 
cant  increase  in  the  age  at  exposure  and  a  decrease  in  7(f)  from  0.038/year  to  0.01 9/year,  R0  changed  only  from  3.3 
to  3.2.  Significant  heterogeneity  was  observed  across  subdistricts  and  schools,  with  R0  ranging  between  1 .7  and 
6.8.  These  findings  are  consistent  with  the  idea  that  the  observed  age  shift  might  be  a  consequence  of  the  demo¬ 
graphic  transition  in  Thailand.  Changes  in  critical  vaccination  fractions,  estimated  by  using  R0,  have  not  accompa¬ 
nied  the  increase  in  age  at  exposure.  These  results  have  implications  for  dengue  control  interventions  because 
multiple  countries  in  Southeast  Asia  are  undergoing  similar  demographic  transitions.  It  is  likely  that  dengue  will 
never  again  be  a  disease  exclusively  of  children. 

basic  reproductive  number;  critical  vaccination  fraction;  demography;  dengue;  force  of  infection;  transmission 
dynamics 


Abbreviation:  Cl,  confidence  interval. 


Dengue  disease  has  been  a  major  public  health  problem  in 
Thailand  over  the  past  50  years.  Dengue  fever  and  dengue 
hemorrhagic  fever  have  traditionally  affected  children,  with 
rare  cases  in  adults,  but  over  the  past  decades,  a  shift  toward 
older  age  groups  has  been  observed  (1-3).  In  the  1980s,  the 
mean  age  of  clinical  cases  was  8-10  years,  but  among  cases 
reported  during  2002-2008,  it  had  increased  to  18  years.  This 
shift  has  occurred  despite  a  roughly  constant  number  of  re¬ 
ported  cases  over  the  same  period.  Similar  shifts  have  been 
reported  in  other  Southeast  Asian  countries  including  Viet¬ 
nam,  Indonesia,  and  Sri  Lanka  (1-5). 

The  following  2  main  hypotheses  have  been  proposed  to 
explain  this  shift:  1 )  Thailand  has  undergone  dramatic  urban¬ 
ization  and  development  over  the  last  decades  that  might  have 


led  to  a  reduction  in  the  transmission  potential  of  dengue 
virus.  Thus,  the  increase  in  the  mean  age  of  cases  could  re¬ 
flect  exposure  later  in  life,  resulting  from  vector  control  pro¬ 
grams  or  decreased  contact  between  humans  and  vectors  (e.g., 
improved  housing  or  increased  air  conditioning  use)  (6-9); 
2)  the  shift  might  be  a  consequence  of  the  demographic  tran¬ 
sition  taking  place  in  Thailand  and  other  Southeast  Asian 
countries.  Mathematical  models  of  disease  transmission  and 
computational  simulations  predict  that  this  transition,  charac¬ 
terized  by  decreasing  birth  rates  and  increasing  life  expect¬ 
ancy,  will  lead  to  an  increase  in  the  average  age  at  infection, 
even  if  the  transmission  potential  of  a  disease  stays  the  same. 
Alternative  hypotheses  include  changes  in  surveillance  pat¬ 
terns  and  changes  in  the  virus  itself  that  might  result  in  more 


353 


Am  J  Epidemiol.  2014;1 79(3):353-360 


Downloaded  from  http://aje.oxfordjoumals.org/  at  AMRMC  on  January  30,  2014 


Report  Documentation  Page 

Form  Approved 

OMB  No.  0704-0188 

Public  reporting  burden  for  the  collection  of  information  is  estimated  to  average  1  hour  per  response,  including  the  time  for  reviewing  instructions,  searching  existing  data  sources,  gathering  and 
maintaining  the  data  needed,  and  completing  and  reviewing  the  collection  of  information.  Send  comments  regarding  this  burden  estimate  or  any  other  aspect  of  this  collection  of  information, 
including  suggestions  for  reducing  this  burden,  to  Washington  Headquarters  Services,  Directorate  for  Information  Operations  and  Reports,  1215  Jefferson  Davis  Highway,  Suite  1204,  Arlington 

VA  22202-4302.  Respondents  should  be  aware  that  notwithstanding  any  other  provision  of  law,  no  person  shall  be  subject  to  a  penalty  for  failing  to  comply  with  a  collection  of  information  if  it 
does  not  display  a  currently  valid  OMB  control  number. 

1 .  REPORT  DATE  2.  REPORT  TYPE 

01  FEB  2014  N/A 

3.  DATES  COVERED 

4.  TITLE  AND  SUBTITLE 

Revisiting  rayong:  shifting  seroprofiles  of  dengue  in  Thailand  and  their 
implications  for  transmission  and  control 

5a.  CONTRACT  NUMBER 

5b.  GRANT  NUMBER 

5c.  PROGRAM  ELEMENT  NUMBER 

6.  AUTHOR(S) 

Rodriguez-Barraquer  I.,  Buathong  R.,  Iamsirithaworn  S.,  Nisalak  A., 
Lessler  J.,  Jarman  R.  G.,  Gibbons  R.  V.,&  Cummings  D.  A. 

5d.  PROJECT  NUMBER 

5e.  TASK  NUMBER 

5f.  WORK  UNIT  NUMBER 

7.  PERFORMING  ORGANIZATION  NAME(S)  AND  ADDRESS(ES) 

US  Army  Institute  of  Surgical  Research,  JBSA  Fort  Sam  Houston,  Texas 

8.  PERFORMING  ORGANIZATION 

REPORT  NUMBER 

9.  SPONSORING/MONITORING  AGENCY  NAME(S)  AND  ADDRESS(ES) 

10.  SPONSOR/MONITOR'S  ACRONYM(S) 

11.  SPONSOR/MONITOR'S  REPORT 
NUMBER(S) 

12.  DISTRIBUTION/AVAILABILITY  STATEMENT 

Approved  for  public  release,  distribution  unlimited 

13.  SUPPLEMENTARY  NOTES 

14.  ABSTRACT 

15.  SUBJECT  TERMS 

16.  SECURITY  CLASSIFICATION  OF:  17.  LIMITATION  OF 

18.  NUMBER  19a.  NAME  OF 

a.  REPORT  b.  ABSTRACT  c.  THIS  PAGE  |J|J 

unclassified  unclassified  unclassified 

8 

Standard  Form  298  (Rev.  8-98) 

Prescribed  by  ANSI  Std  Z39-18 


354  Rodriguez-Barraquer  et  al. 


clinical  disease  among  adults  or  reinfection  of  those  immune 
to  previously  circulating  viruses. 

For  multiserotype  diseases  like  dengue,  for  which  the  fac¬ 
tors  associated  with  symptomatic  disease  have  not  been  fully 
characterized,  it  is  not  clear  that  the  age  distribution  of  clin¬ 
ically  apparent  cases  represents  the  age  distribution  of  in¬ 
fection.  Secondary  infection  by  a  heterotypic  serotype  has 
consistently  been  found  to  be  a  main  risk  factor  for  develop¬ 
ing  severe  disease  after  infection  (10),  but  the  associations 
among  age,  risk  of  infection,  and  clinical  disease  are  not 
well  understood.  In  such  instances,  inference  from  the  age 
distribution  of  cases  can  be  misleading,  and  age  at  exposure 
can  be  determined  only  through  age-stratified  serosurveys. 

In  December  2010,  we  conducted  an  age-stratified,  school- 
based  serological  survey  in  Mueang  Rayong  district,  Thailand. 
This  study  was  designed  to  replicate  a  serosurvey  conducted 
in  the  same  location  in  1980/1981  to  the  extent  possible  (11, 
12).  By  using  data  from  the  2  studies,  we  estimated  serotype- 
specific  forces  of  infection  (k)  and  basic  reproductive  num¬ 
bers  ( R0 ,  a  measure  of  transmission  potential)  for  each  of  the 
dengue  serotypes  for  the  periods  1969-1980  and  1993-2010. 
Our  results  suggest  that,  despite  a  50%  decrease  in  the  force 
of  infection,  R0  (and  hence  transmission  potential)  has  not 
changed  significantly.  These  results  suggest  that  factors  other 
than  a  decrease  in  the  transmission  potential  might  be  driving 
the  epidemiologic  shift  of  dengue  in  Thailand. 

MATERIALS  AND  METHODS 

Rayong  Province  is  located  in  southeast  Thailand  and  is 
hyperendemic  for  dengue.  It  has  historically  reported  1  of 
the  highest  rates  of  dengue  in  Thailand  (13).  The  reported 
dengue  incidence  was  183.6  per  100,000  population  (116,947 
cases)  in  2010  and,  of  these,  52.0%  (60,770)  were  cases  of 
dengue  hemorrhagic  fever  (14).  We  used  data  from  age- 
stratified  serosurveys  conducted  30  years  apart  in  the  Mueang 
Rayong  district. 

1980/1981  Surveys 

In  1980  and  1981,  2  school-based  age- stratified  serosurveys 
were  conducted  in  the  municipal  area  of  Mueang  Rayong, 
Thailand,  as  part  of  a  then  ongoing  cohort  study.  Details  about 
this  study  have  been  published  elsewhere  (11,  12).  Dengue 
serostatus  was  determined  by  using  a  single  dilution  (1:30) 
neutralization  test.  Past  exposure  to  dengue  virus  (to  1  or 
more  serotypes)  was  assumed  if  there  was  70%  or  more  pla¬ 
que  reduction  to  1  or  more  dengue  serotypes  at  a  1:30  dilu¬ 
tion.  In  a  subsample,  the  endpoint  titer  needed  to  achieve 
50%  plaque  reduction  was  calculated  from  2  or  more  dilu¬ 
tions.  Data  from  the  1980  serosurvey  were  previously  ana¬ 
lyzed  by  Ferguson  et  al.  (15)  to  produce  estimates  of  k  and  R0. 

201 0  Survey 

In  December  2010,  we  conducted  a  new  school-based, 
age-stratified  serosurvey  in  Mueang  Rayong  district.  This 
serosurvey  was  not  limited  to  the  municipal  area  but  in¬ 
cluded  all  of  the  subdistricts  within  Mueang  Rayong.  A  prob¬ 
abilistic  sample  of  25  schools  and  classes  were  selected  from 


administrative  data  provided  by  the  Rayong  school  district 
administration,  and  children  aged  4-18  years  were  sampled. 
Dengue  serostatus  was  determined  by  using  the  single  dilu¬ 
tion  neutralization  test  against  the  4  dengue  serotypes  at  a  sin¬ 
gle  1:10  dilution.  The  decision  to  use  a  single  1:10  dilution 
instead  of  a  1:30  dilution  was  based  on  improved  sensitivity 
when  characterizing  historical  infections  (R.J.,  unpublished 
data,  2009)  and  a  lack  of  clarity  in  the  original  survey  as  to 
the  precise  methods  used  for  each  reported  sample.  Past  ex¬ 
posure  to  dengue  virus  was  defined  as  70%  or  more  plaque 
reduction  to  1  or  more  dengue  serotypes.  We  used  a  higher 
(90%)  plaque  reduction  cutoff  to  differentiate  children  with 
monotypic  versus  multitypic  immunity  in  order  to  optimize 
the  specificity  of  the  assay  (R.J.,  unpublished  data,  2009). 
We  also  assessed  several  additional  plaque  reduction  cut- 
points. 

By  using  the  1980/1981  and  2010  data  sets,  we  estimated 
the  average  and  serotype-specific  forces  of  infection  kit)  and 
basic  reproductive  number  R0  in  the  Mueang  Rayong  municipal 
area.  To  characterize  spatial  heterogeneity  in  transmission,  we 
also  used  the  data  from  the  2010  survey  to  estimate  subdistrict- 
and  school-specific  kit)  and  Rq.  All  statistical  analyses  were 
performed  in  R  program,  version  2.14.0  (R  Foundation  for 
Statistical  Computing,  Vienna,  Austria)  and  JAGS,  version 
3.1.0  (16). 

Estimating  the  force  of  infection 

The  force  of  infection  (X)  is  used  to  characterize  transmis¬ 
sion  hazard  in  a  given  setting  and  is  defined  as  the  rate  at 
which  susceptible  individuals  are  infected.  Age- stratified 
serosurveys  can  be  used  to  estimate  X.  By  assuming  that  the 
risk  of  infection  does  not  vary  with  age  and  that  immunity  is 
lifelong,  the  difference  in  seroprevalence  between  subjects  a 
and  fl+l  year  of  age  can  be  attributed  to  transmission  that  oc¬ 
curred  between  a  and  a  +  1  year  ago.  Methods  to  estimate 
forces  of  infection  from  seroprevalence  data  have  been  de¬ 
scribed  elsewhere  (15,  17-20). 

Given  the  age  ranges  in  the  seroprevalence  surveys  avail¬ 
able,  we  were  able  to  estimate  X(t)  for  the  periods  1969-1980 
(from  the  1980/1981  serosurvey)  and  1993-2010  (from  the 
2010  serosurvey).  We  estimated  forces  of  infection  averaged 
across  all  serotypes  (A,(f)),  as  well  as  for  each  serotype  (X, ■(?))■ 
We  fit  models  assuming  a  constant  force  of  infection  through¬ 
out  time  and  models  that  allowed  for  yearly  time- varying  forces 
of  infection  (piecewise  constant).  To  produce  subdistrict-  and 
school-specific  estimates  of  X(t),  we  also  fit  a  random-effects 
model. 

Estimating  the  basic  reproductive  number  ( R0 ) 

R0  is  the  number  of  secondary  infections  generated  by  a  pri¬ 
mary  infection  in  a  completely  susceptible  population  during 
the  entire  duration  of  infectiousness  and,  as  such,  is  a  measure 
of  transmission  potential.  For  vector-borne  diseases,  multiple 
R0  can  be  estimated  (e.g.,  human-to-mosquito,  mosquito-to- 
human,  human-to-human).  Flere,  we  estimate  the  number  of 
secondary  human  infections  that  arise  from  1  infectious  hu¬ 
man  as  mediated  through  the  vector,  as  proposed  by  Ferguson 
et  al  in  1999  (15).  As  with  kit),  we  provide  average  estimates 
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of  R0  across  all  serotypes,  as  well  as  serotype-specific  esti¬ 
mates.  Our  estimates  of  R0  assume  that  individuals  can  be  in¬ 
fected  by  up  to  2  serotypes  and  are  afterward  immune  to  further 
infection. 

At  equilibrium,  R0  is  related  to  the  critical  vaccination  frac¬ 
tion  (  V)  by  the  equation  V  =  1  —  1  /R0  and  hence  gives  insight 
into  the  level  of  control  that  is  required  to  eventually  block 
transmission  (21).  We  estimated  R0  for  Mueang  Rayong  by 
using  our  Ml)  estimates  and  data  from  the  1980  and  2010 
Thailand  censuses  (22,  23). 

Full  details  on  all  methods  can  be  found  in  Web  Appendix  1 
and  Web  Figure  1,  available  at  http://aje.oxfordjournals.org/. 
The  protocol  for  this  study  was  reviewed  and  approved  by 
the  institutional  review  boards  of  the  Thai  Ministry  of  Public 
Health  (Nonthaburi,  Thailand)  and  the  Johns  Hopkins 
Bloomberg  School  of  Public  Health  (Baltimore,  Maryland). 

RESULTS 

The  1980  and  1981  surveys  contained  data  from  1,009  and 
1,235  children  aged  0-10  years  and  0-7  years,  respectively. 
The  2010  serosurvey  included  data  from  1,736  children  aged 


4-18  years,  of  whom  898  lived  in  the  Mueang  Rayong  mu¬ 
nicipal  area  (area  of  the  1980  and  1981  studies).  We  excluded 
from  analysis  89  children  aged  4  or  5  years  from  the  2010 
data  set  because  of  their  low  participation  rate  (20%-30%) 
compared  with  all  other  age  groups  (>50%). 

Overall,  59%  of  the  sampled  children  showed  evidence  of 
prior  exposure  to  dengue  virus  in  the  2010  study.  Figures  1 
and  2  show  the  age-specific  seroprevalence  in  the  2010  sero¬ 
survey  compared  with  the  1980  serosurvey.  According  to  the 
1980  survey,  96%  of  the  population  had  been  exposed  to 
dengue  by  age  11  years.  The  2010  survey  estimated  that 
only  65%  (95%  confidence  interval  (Cl):  61%,  75%)  of  chil¬ 
dren  aged  1 1  years  had  been  exposed  and  that,  at  age  18  years, 
18%  (95%  Cl:  10%,  26%)  of  the  population  remained  sus¬ 
ceptible  to  dengue  virus. 

Force  of  infection 

The  estimated  time-constant  force  of  infection  (A,(f))  for  the 
period  1969-1980  was  0.038  per  year  (95%  Cl:  0.036,  0.039). 
This  means  that,  on  average,  3.7%  (1  —  e-0'038)  of  suscepti¬ 
ble  individuals  were  infected  by  each  serotype  each  year 
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Figure  1 .  Age-specific  dengue  seroprevalence  in  Mueang  Rayong  district,  Thailand,  according  to  the  A)  1 980,  and  B)  201 0  serosurveys.  Numbers 
on  the  top  of  each  bar  indicate  the  number  of  samples  in  each  age  group.  Whereas  the  1 980  serosurvey  showed  that  96%  of  the  population  had  been 
exposed  to  dengue  by  age  1 1  years,  the  2010  serosurvey  showed  that  only  82%  of  the  population  was  exposed  by  age  18  years.  DENV-1,  DENV-2, 
DENV-3,  and  DENV-4  refer  to  dengue  virus  serotypes  1-4,  respectively,  and  >2  DENV  refers  to  2  or  more  dengue  virus  serotypes. 
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Figure  2.  Age-specific  dengue  seroprevalence  in  Mueang  Rayong 
district,  Thailand,  according  to  the  1980  and  2010  serosurveys.  Solid 
lines  show  the  fits  of  the  constant  X(f)  (serotype-specific  forces  of  infec¬ 
tion)  models.  Dashed  lines  show  the  fits  of  the  time-varying  X(t)  model. 
Bars  represent  95%  confidence  intervals.  Results  show  a  significant 
increase  in  the  age  at  exposure  to  dengue  virus.  DENV-1 ,  DENV-2, 
DENV-3,  and  DENV-4  refer  to  dengue  virus  serotypes  1-4,  respec¬ 
tively,  and  >2  DENV  refers  to  2  or  more  dengue  virus  serotypes. 


Appendix  2,  Web  Tables  1  and  2,  and  Web  Figure  2)  also 
suggest  decreases  in  X(t)  that  far  exceed  changes  in  Rq. 

Serotype-specific  estimates  of  X  and  R0 

Figure  3B  shows  our  serotype-specific  estimates  of  the 
basic  reproductive  number  and  force  of  infection  for  the  2  time 
periods.  These  estimates  are  also  consistent  with  greater  reduc¬ 
tions  in  the  force  of  infection  across  all  serotypes  (Figure  3  A) 
than  in  the  basic  reproductive  number.  Interestingly,  dengue 
virus  serotype  2  seems  to  continue  to  cause  the  highest  hazard 
of  infection  in  Mueang  Rayong  district.  Between  1993  and 
2010,  Rq  and  the  critical  vaccination  fraction  for  dengue  virus 
serotype  2  were  3.7  (95%  Cl:  3.5,  3.8)  and  73%,  respectively, 
which  was  significantly  higher  than  the  corresponding  values 
for  the  other  serotypes. 

As  expected,  our  serotype-specific  estimates  of  kit)  for  the 
period  1969-1980  are  similar  to  those  produced  by  Ferguson 
et  al.  (15)  for  the  period  1969-1979  (using  the  1980  data  set 
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during  this  period.  In  contrast,  the  estimated  time-constant 
X(f)  for  the  period  1993-2010  was  0.019  per  year  (95%  Cl: 
0.017,  0.020)  for  the  municipal  area,  ranging  between  0.007 
and  0.025  among  schools  in  the  municipality.  These  results 
suggest  a  50%  decrease  in  X  between  the  2  periods  (Figure  2). 

Allowing  for  time-varying  X(t)  (piecewise  constant)  (Figure  2) 
significantly  improved  the  fit  of  the  model  to  the  1980/1981  age- 
specific  seroprevalence  data  (likelihood  ratio  test,  P  <  0.001) 
but  did  not  significantly  improve  the  fit  to  the  2010  data  (like¬ 
lihood  ratio  test,  P  =  0.1).  Results  of  this  time-varying  model 
suggest  that  transmission  intensity  was  higher  in  1969  and 
1970  than  in  the  rest  of  the  estimated  period  but  should  be 
interpreted  with  caution  because  of  the  cross-sectional  nature 
of  the  data. 


Basic  reproductive  number  and  critical  vaccination 
fraction 

By  using  the  estimates  of  X(t)  and  age-specific  census  data 
from  the  1980  and  2010  censuses,  we  estimated  the  average 
Rq  to  have  been  3.3  (95%  Cl:  3.1,  3.4)  for  the  period  1 969— 
1980,  corresponding  to  a  critical  vaccination  fraction  of  70%. 
The  estimated  R0  for  the  period  1993-2010  was  3.2  (95%  Cl: 
2.7,  3.3)  for  the  municipal  area,  but  ranged  from  1.8  to  4.3 
among  schools  in  this  area.  These  R0  estimates  correspond 
to  a  critical  vaccination  fraction  of  69%  (range,  44%-77%). 
Thus,  our  results  suggest  that,  despite  a  50%  decrease  in  k(t) 
over  the  last  30  years,  Rq  and  the  critical  vaccination  fraction 
have  not  changed.  Results  from  sensitivity  analyses  (Web 
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Figure  3.  Serotype-specific  estimates  of  A)  X(f)  (serotype-specific 
forces  of  infection),  and  B)  R0  (basic  reproductive  number)  for  dengue 
for  the  periods  1969-1980  and  1993-2010  in  Mueang  Rayong  district, 
Thailand.  DENV-1,  DENV-2,  DENV-3,  and  DENV-4  refer  to  dengue 
virus  serotypes  1-4,  respectively.  Although  the  force  of  infection  has 
decreased  by  approximately  50%,  changes  in  R0  have  been  signifi¬ 
cantly  smaller. 
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only).  Our  RQ  estimates  are  slightly  higher  because  of  the  as¬ 
sumptions  we  made  regarding  interaction  between  serotypes 
(Web  Table  1). 

Spatial  heterogeneity  in  X  and  R0 

Figure  4  shows  our  subdistrict-specific  estimates  of  X(/j 
for  the  period  1993-2010.  The  sample  included  schools  in 
9  of  the  15  subdistricts  in  Mueang  Rayong;  therefore,  we 
were  able  to  estimate  X(f)  only  for  these  9  locations. 

Our  results  suggest  that  there  is  significant  heterogeneity  in 
transmission  intensity  within  the  Mueang  Rayong  district. 
The  between-subdistrict  variance  in  X(/j  was  estimated 
to  be  0.50  (95%  Cl:  0.27,  0.54).  Relative  to  Tha  Pradu,  the 
subdistrict  that  includes  the  central  urban  area  of  Mueang 
Rayong,  the  X(f)s  of  subdistricts  Choeng  Noen  and  Map  Ta 
Phut  were  significantly  higher  during  the  period  1993-2010 
(relative  X(f)s  of  1.35  (95%  Cl:  1.12,  1.63)  and  1.91  (95%  Cl: 
1.55,  2.30),  respectively),  and  X(/j  in  Ka  Chet  was  signifi¬ 
cantly  lower  (relative  X(Y)  of  0.49  (95%  Cl:  0.34,  0.67)). 
Estimates  of  subdistrict-specific  X(/js  and  relative  X(/js  are 
provided  in  Web  Table  3. 

Heterogeneity  in  transmission  was  also  significant  when 
looking  at  finer  spatial  scales.  The  between-school  variance 
in  X(/j  was  estimated  to  be  0.42  (95%  Cl:  0.20,  0.79)  for 
schools  located  in  the  municipality,  and  X(f)  ranged  from 


Figure  4.  Estimated  X(f)  (serotype-specific  forces  of  infection)  for  the 
9  subdistricts  of  Mueang  Rayong,  Thailand,  represented  in  the  2010 
serological  study.  Locator  map  on  the  upper  left  corner  shows  the  lo¬ 
cation  of  Mueang  Rayong  district  (black  point)  within  Thailand.  Results 
show  significant  heterogeneity  in  X(f)  among  subdistricts.  A)  Tha 
Pradu,  B)  Choeng  Noen,  C)  Noen  Phra,  D)  Taphong,  E)  Phe, 
F)  Klaeng,  G)  Na  Ta  Khuan,  H)  Ka  Chet,  and  I)  Map  Ta  Phut.  Credible 
intervals  of  our  estimates  are  provided  in  Web  Table  1 ,  available  at  http:// 
aje  .oxfordjou  rnals.org/. 


0.009  to  0.024  among  these  schools.  Estimates  of  X(/j  ranged 
between  0.007  and  0.044  when  all  schools  in  the  district  were 
considered. 

As  expected  from  the  heterogeneity  in  kit),  R0  and  V 
were  also  heterogeneous  across  subdistricts  (Figure  5). The 
mean  R0  in  Mueang  Rayong  district  estimated  by  using  this 
model  was  3.2  (95%  Cl:  2.52,  4.04)  but  ranged  from  1.82 
(95%  Cl:  1.52,  2.20)  in  Ka  Chet  to  5.13  (95%  Cl:  4.38, 
5.96)  in  Map  Ta  Phut.  This  implies  that  the  mean  critical  vac¬ 
cination  fraction  can  range  from  45%  to  81%  in  different  sub¬ 
districts.  School-specific  estimates  of  R0  ranged  between  1.71 
(V=0.41)  and  6.76  (F=0.85)  and  are  presented  in  Web 
Figure  3. 

Though  it  is  likely  that  part  of  this  observed  heterogeneity 
could  be  explained  by  demographic  or  socioeconomic  fac¬ 
tors,  none  of  the  predictors  measured  in  this  study  (sex, 
number  of  household  members,  access  to  electricity,  vehicle 
ownership,  availability  of  running  water,  residence  time  at 
address)  showed  a  significant  association  when  included  in 
the  random-effects  model.  Interestingly,  Map  Ta  Phut  is  a 
major  industrial  subdistrict  in  the  region  with  the  highest 
rates  of  immigration  and  mobility  of  the  population. 


Figure  5.  Heterogeneity  of  R0  (basic  reproductive  number)  and  crit¬ 
ical  vaccination  fraction  estimates  between  subdistricts  of  Mueang 
Rayong,  Thailand.  Estimates  obtained  from  the  2010  serological 
study.  Solid  line  represents  the  equation  V  (critical  vaccination  frac¬ 
tion)  =  1-1  /R0.  Squares  represent  the  I/and  R0  estimates  for  the  9 
subdistricts  included  in  the  serosurvey.  A)  Tha  Pradu,  B)  Choeng 
Noen,  C)  Noen  Phra,  D)  Taphong,  E)  Phe,  F)  Klaeng,  G)  Na  Ta  Khuan, 
H)  Ka  Chet,  and  I)  Map  Ta  Phut.  Darker  shades  represent  higher  R0. 
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DISCUSSION 

Thailand  and  several  other  Southeast  Asian  countries  have 
experienced  a  dramatic  shift  in  the  age  distribution  of  dengue 
cases  over  the  last  30  years.  Traditionally,  severe  dengue  af¬ 
fected  young  children,  but  during  recent  years,  cases  have  be¬ 
come  increasingly  common  among  adolescents  and  young 
adults  (3,  5)  despite  a  roughly  constant  incidence.  Our  results 
show  that  underlying  this  shift  there  have  been  significant 
changes  in  the  age  of  those  exposed  and  patterns  of  exposure 
to  dengue  virus. 

Although  several  studies  have  reported  shifts  in  the  age  of 
dengue  cases,  all  of  these  have  been  based  on  case  data  alone 
and  can  therefore  be  biased  because  of  age  dependence  of 
clinical  disease  (24,  25).  Our  results,  from  age-stratified  sero¬ 
logical  studies  conducted  in  the  same  location  30  years  apart, 
suggest  a  decrease  in  the  force  of  infection  of  approximately 
50%,  with  a  much  smaller,  nonsignificant  change  in  R0. 
These  findings  are  consistent  with  case  data-based  results 
reported  previously  by  several  coauthors  (3).  By  using  age- 
specific  incidence  data,  Cummings  et  al.  (3)  estimated  k(t) 
and  R0  for  the  period  1980-2005  for  72  provinces  and  showed 
a  significant  decrease  in  k(f)  accompanied  by  a  slighter  change 
in  R0.  These  findings  imply  that,  despite  urbanization  and  de¬ 
velopment  in  Mueang  Rayong  district  during  the  last  30  years, 
the  transmission  potential  and  critical  vaccination  fraction  of 
dengue  have  not  significantly  changed. 

Our  observations  are  consistent  with  the  shift  in  age  of  ex¬ 
posure  to  dengue  being  driven  by  changes  in  the  age  structure 
of  the  Thai  population  over  the  last  30  years  (Web  Figure  4) 
(3).  The  idea  that  changes  in  the  age  structure  of  populations 
can  lead  to  changes  in  the  transmission  hazard  and  age  distri¬ 
bution  of  cases  has  been  discussed  for  several  infectious  dis¬ 
eases  (17,  26).  Similarly,  the  impact  of  population  growth  on 
R0  has  also  been  discussed  (27,  28).  Decreases  in  birth  rate 
and  increases  in  life  expectancy  lead  to  a  decrease  in  both 
the  rate  at  which  susceptible  individuals  enter  the  population 
and  the  rate  at  which  immune  individuals  are  removed  from 
it.  This  results  in  smaller  households,  an  increased  probabil¬ 
ity  of  contact  with  older  (immune)  people,  and,  in  the  case  of 
vector-borne  diseases,  an  increase  in  the  likelihood  that  infec¬ 
tious  bites  will  be  delivered  to  an  already  immune  person.  Be¬ 
cause  of  this  reduced  hazard,  more  individuals  remain 
susceptible  to  primary  or  secondary  infection  into  adulthood, 
and  infections  come  from  a  broader  range  of  age  classes. 

Although  our  results  suggest  that  the  demographic  transi¬ 
tion  might  be  a  major  driving  force  behind  the  observed 
changes  in  the  epidemiology  of  dengue  in  Mueang  Rayong 
district,  they  do  not  preclude  that  other  factors  might  also 
play  a  role.  Mueang  Rayong  has  experienced  major  changes 
over  the  last  30  years,  becoming  a  tourist  and  industrial  center 
in  central/eastern  Thailand.  Mueang  Rayong  municipality 
has  undergone  progressive  urbanization;  its  population  has 
grown  from  38,000  in  1980  to  58,000  in  2010,  and  its  density 
has  increased  from  2,190  people/km2  to  more  than  3,400  peo¬ 
ple/km2  (29).  The  mean  size  of  households  has  decreased 
from  5.6  in  1980  to  3.1  in  2000,  and  access  to  electricity 
and  tap  water  has  become  almost  universal  (30).  Further¬ 
more,  Thailand  reinforced  its  vector  control  programs  and  ed¬ 
ucational  campaigns  in  the  1990s  (31,  32).  Though  each  of 


these  factors  may  contribute  to  changes  in  the  serological  pat¬ 
tern  of  dengue  infection,  major  impacts  are  inconsistent  with 
the  small  change  in  R0  over  these  last  30  years.  Extensive  lon¬ 
gitudinal  data  (serological,  demographic,  and  entomological) 
would  be  necessary  to  fully  characterize  the  relative  contribu¬ 
tion  of  entomological  versus  demographic  factors  to  the  cur¬ 
rent  epidemiologic  shift  observed  in  dengue  in  Thailand  and 
other  Southeast  Asian  countries. 

We  present  the  first  estimates  in  30  years  of  the  basic  repro¬ 
ductive  number  of  dengue  derived  from  serological  data.  Our 
estimates  suggest  that  R0  of  dengue  virus  serotype  2  is  signif¬ 
icantly  higher  than  that  of  other  serotypes,  implying  that  it 
would  require  a  higher  level  of  vaccination  to  control  trans¬ 
mission.  However,  though  our  findings  show  some  degree  of 
heterogeneity  in  transmission  potential  among  serotypes, 
other  sources  of  heterogeneity  seem  to  be  more  important. 
In  particular,  our  results  suggest  that  the  critical  vaccination 
threshold  might  vary  by  as  much  as  45%  in  different  loca¬ 
tions  within  the  same  district,  and  this  difference  might  be 
greater  when  looking  at  finer  spatial  scales.  Spatial  heteroge¬ 
neity  in  transmission  potential  has  been  previously  described 
for  dengue  and  other  vector-borne  diseases  (33-35)  and  im¬ 
plies  that  higher  vaccination  coverages  or  targeted  vaccina¬ 
tion  might  be  required  to  effectively  control  transmission. 
Further  characterization  of  the  extent  of  heterogeneity  will 
be  fundamental  to  designing  proper  vaccination  strategies. 

Although  the  plaque  reduction  neutralization  assay  is  con¬ 
sidered  the  “gold  standard”  to  characterize  immunity  after 
dengue  exposure,  optimal  cutpoints  to  classify  immune  re¬ 
sponses  for  seroprevalence  studies  have  not  been  established. 
Optimal  neutralization  cutpoints  for  the  single  dilution  neu¬ 
tralization  test  used  in  this  study  are  also  lacking,  particularly 
for  differentiating  monotypic  versus  multitypic  immunity. 
We  present  results  from  70%  and  90%  neutralization  based 
on  data  (R.J.,  unpublished  data,  2009)  suggesting  that  these 
might  be  the  optimal  cutpoints  in  terms  of  sensitivity  and  spec¬ 
ificity,  but  we  performed  sensitivity  analyses  for  different 
neutralization  levels  (Web  Figure  5). 

Although  estimates  of  X(t)  from  age-stratified  serological 
data  are  robust  to  changes  in  demography,  this  method  of  es¬ 
timating  R0  does  make  the  simplifying  assumption  that  the 
population  structures  during  each  period  (i.e.,  1969-1980 
and  1993-2010)  were  stable  and  equal  to  the  population 
structures  in  1980  and  2010,  respectively.  Properly  account¬ 
ing  for  changes  in  demography  would  require  incorporating 
detailed  longitudinal  population  data  (mortality,  births,  mi¬ 
gration)  into  the  model  and  is  beyond  the  scope  of  our  study. 
Sensitivity  analyses  in  which  we  used  multiple  population 
structures  for  each  time  period  show  that  our  findings  are  ro¬ 
bust  to  this  assumption,  because  the  estimated  reduction  in 
/,( 1)  between  the  2  periods  was  always  much  larger  than  the 
reduction  in  R0  (Web  Table  1). 

Estimates  of  k(t)  and  R0  from  seroprevalence  data  are  also 
sensitive  to  model  choice  and  to  assumptions  made  regarding 
mixing  of  the  population  and  interactions  between  serotypes, 
among  others.  Estimates  of  R0  are  also  sensitive  to  assump¬ 
tions  made  regarding  force  of  infection  during  unobserved 
periods  (e.g.,  before  1969  and  during  1981-1993).  Our  esti¬ 
mates  assume  homogeneous  mixing  and  full  immunity  after 
2  infections  by  heterologous  serotypes  and  ignore  short-lived 
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cross-protection  between  serotypes.  They  also  assume  a 
constant  hazard  of  infection  prior  to  1969  (equal  to  the  con¬ 
stant  hazard  in  the  period  1969-1980).  Results  from  sensitiv¬ 
ity  analyses  suggest  that  our  findings  are  also  robust  to  these 
assumptions  (Web  Tables  1  and  2  and  Web  Figure  2). 

These  findings  have  important  implications  in  the  design 
and  implementation  of  interventions  and  highlight  the  rele¬ 
vance  of  considering  changes  in  demography  as  drivers  of 
changes  in  age  distribution  of  infectious  diseases.  If  the  de¬ 
mographic  transition  is  in  fact  a  major  driving  force  behind 
the  observed  changes  in  the  epidemiology  of  dengue  in  Thai¬ 
land  and  Southeast  Asia,  it  is  very  likely  that  dengue  will 
never  again  be  a  disease  exclusively  of  children.  Epidemio¬ 
logic  surveillance,  control  programs,  and  management  of  den¬ 
gue  will  need  to  expand  and  focus  on  the  detection,  treatment, 
and  prevention  of  both  adult  and  childhood  cases. 
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